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LONG-TERM  GOALS 

The  long-term  goal  of  this  project  is  to  fully  utilize  microwave  backscattering  data  sets  that  have  been 
collected  on  the  ocean  over  the  last  two  decades.  The  prime  objective  is  to  better  understand  the 
formation  of  breaking  waves  on  the  ocean  and  their  effect  on  microwave  backscatter. 

SCIENTIFIC  OBJECTIVES 

The  scientific  objectives  of  this  research  are  to  investigate  experimentally  wave  shadowing  and 
modulation  in  low-grazing-angle  backscatter  from  the  ocean. 

APPROACH 

Our  approach  is  to  reanalyze  data  sets  taken  with  Doppler  radars  over  the  last  two  decades  in  an 
attempt  to  extract  additional  information  on  microwave  backscatter  from  the  ocean.  These  data  sets 
have  been  taken  on  a  variety  of  platforms  from  airplanes  and  towers  to  blimps  and  ships  using  a  variety 
of  coherent  microwave  radars.  In  the  particular  work  that  has  been  carried  out  in  the  past  year,  we 
have  looked  for  evidence  of  shadowing  of  microwave  backscatter  by  ocean  waves  and  have  applied 
bound/breaking  wave  theory  to  observed  phase  differences  between  modulated  received  power  and 
scatterer  velocities. 

WORK  COMPLETED 

Several  data  sets  collected  over  the  past  two  decades  were  examined  to  see  if  they  supported  a 
geometric  shadowing  picture  of  low- grazing-angle  backscatter  from  the  ocean.  We  may  consider 
effects  of  the  possible  types  of  rough-surface  shadowing  on  the  modulation  of  backscattered  power  by 
ocean  waves  as  shown  in  Figure  1.  Incoming  rays  are  indicated  by  slanting  lines  from  the  left.  Figures 
la  and  c  show  geometric  shadowing  at  low  and  high  incident  power  levels.  Increasing  the  incident 
power  can  only  increase  the  return  in  the  illuminated  areas.  By  the  very  definition  of  geometric 
shadowing,  the  change  in  power  level  cannot  affect  the  regions  that  are  not  illuminated  so  the 
probability  of  observing  signals  at  the  noise  level  cannot  change.  Also  note  that  geometric  shadowing 


1 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

30  SEP  2011  2' REPORT  TYPE 

3.  DATES  COVERED 

00-00-2011  to  00-00-2011 

4.  TITLE  AND  SUBTITLE 

Analysis  of  Coherent  Microwave  Data  Collected  on  the  Ocean  over  Two 
Decades 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Washington, Applied  Physics  Laboratory, 1013  NE  40th 

St, Seattle, WA, 98105-6698 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

ARSTRATT 

1 8 .  NUMBER  1 9a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  aS 

unclassified  unclassified  unclassified  Report  (SAR) 

8 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


is  a  ray  theory  so  the  polarization  of  the  incident  radiation  is  irrelevant.  By  contrast,  Figures  lb  and  d 
show  cases  of  partial  shadowing  at  low  and  high  incident  power  levels.  In  this  case,  backscatter  can 
come  from  the  areas  that  are  geometrically  shadowed  but  may  be  reduced  either  because  the  incident 
field  is  reduced  or  because  scatterers  in  this  region  are  weaker  than  in  the  illuminated  areas  due  to  their 
modulation  by  waves.  Both  of  these  processes  may  be  polarization  dependent.  If  the  incident  power 
level  is  low,  the  radar  may  still  observe  signals  at  the  noise  level  in  areas  behind  a  crest  as  shown  in 
Figure  lb.  In  practice,  the  difference  between  the  modulations  of  received  power  diagrammed  in 
Figures  la  and  b  may  be  difficult  to  distinguish.  However,  for  partial  shadowing,  an  increase  in  the 
incident  power  level  may  raise  the  backscatter  from  behind  the  crests  above  the  noise  level  as  shown  in 
Figure  Id.  Thus  either  changing  the  incident  power  level  or  noting  polarization  differences  may  allow 
differentiation  between  geometric  and  partial  shadowing. 


Figure  1.  Diagrams  of  geometric  and  partial  shadowing  at  low  and  high  incident  power  levels. 


Pr  is  received  power,  N  is  the  radar  noise  level  indicated  by  the  dashed  line,  and  r|  is  the  sea  surface,  a) 
geometric  shadowing,  low  incident  power;  b)  partial  shadowing,  low  incident  power;  c)  geometric 
shadowing,  high  incident  power;  d)  partial  shadowing,  high  incident  power. 


For  all  of  the  data  collected  with  our  radars,  we  computed  the  signal-to-noise  ratio,  SNR: 


SNR  = 


p;-  n 

N 


N 


where  Pr'  is  the  received  power  from  the  sea  surface  plus  noise.  We  then  determined  the  fraction  of  the 
data  from  various  runs  that  had  SNR<1.  We  will  call  this  fraction  the  concealed  fraction  since  the 
signal,  if  any,  is  smaller  than  the  system  noise  and  therefore  concealed  from  detection.  Each  sample  of 
received  power  was  averaged  only  over  the  time  that  the  radars  took  to  collect  them,  either  262  msec 
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or  41  msec.  Each  data  set  contained  in  excess  of  300  data  points  at  each  range  bin,  each  of  which 
yielded  a  value  of  SNR.  The  number  of  SNR  values  less  than  one  was  then  determined  and  divided  by 
the  total  number  of  samples  for  that  range  bin  for  that  run. 

Average  SNR  values  in  dB  divided  by  100  and  the  concealed  fractions  at  various  ranges  for  HH 
polarized  data  collected  on  an  airship  in  1995  are  shown  in  Figure  2.  The  antenna  height  was  240  m  so 
the  grazing  angles  for  the  ranges  shown  are  from  25.2°  to  71.6°.  In  this  range  of  grazing  angles, 
shadowing  is  not  expected  and  the  concealed  fraction  becomes  non-zero  only  where  the  incident  power 
is  low.  This  shows  that,  as  expected,  weak  scatterers  can  produce  very  low  backscatter  even  in  the 
absence  of  shadowing  if  the  incident  power  level  is  low. 


Figure  2.  Data  collected  on  the  airship  in  1995  with  HH  polarization  for  grazing  angles  between 
25.2°  and  71.6°.  The  solid  curve  is  the  average  signal-to-noise  ratio,  SNR,  in  dB  divided  by  100  and 
the  asterisks  show  the  fraction  of  SNR  values  having  SNR<1,  the  concealed  fraction.  Wind  12.6  m/s 
from  1 73°T,  antenna  boresite  55°  looking  toward  180°T.  Range  resolution  in  7.5  m. 

Figure  3  shows  VV  polarized  data  taken  on  the  R/V  Thompson  in  2008  at  very  low  grazing  angles. 
Now  the  concealed  fraction  is  significant  even  where  the  incident  power  is  high.  This  could  indicate 
that  some  geometric  shadowing  is  taking  place.  However,  the  concealed  fraction  exhibits  much 
different  behavior  than  predicted  by  Milder  (2004)  whose  result  is  shown  as  the  dashed  line  in  Figure 
3.  This  was  computed  as  the  fraction  shadowed 

where  f  is  Milder’s  fraction  illuminated,  which  he  gives  in  his  Eq.3.7  as 

f(cc)  =  V2ri(a/cr) 

where  a  is  the  tangent  of  the  grazing  angle  and  a  is  the  rms  slope  of  the  rough  surface.  Milder 
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says  that  this  result  is  valid  when  a  «  cr.  We  show  values  when  af  a  <  1/V2m .  For  a  we  have  used 
the  wind- speed  dependent  results  of  Cox  and  Munk  (1954)  for  a  slick-covered  surface,  which  makes  o 
as  small  as  possible  and  therefore  fs  as  small  as  possible.  We  see  that  even  using  this  small  rms  slope, 
the  predicted  fraction  shadowed  is  generally  larger  than  observed  concealed  fraction  except  at  grazing 
angles  in  the  near  range  where  the  theory  says  that  geometrical  shadowing  should  not  occur. 
Furthermore,  in  the  data,  the  concealed  fraction  decreases  with  increasing  wind  speed  while  the  theory 
predicts  the  opposite.  If  the  theory  is  correct,  and  other  theories  of  geometric  shadowing  are  very 
similar,  then  again  geometrical  shadowing  does  not  appear  to  explain  the  observations. 

In  the  data  discussed  so  far,  the  decrease  in  incident  power  for  any  particular  set  of  environmental  and 
radar  parameters  was  always  correlated  with  a  decrease  in  grazing  angle.  Thus  it  could  be  argued  that 
the  increase  in  concealed  fraction  at  far  ranges  could  be  due  to  the  decreasing  grazing  angle,  and 
subsequent  increasing  shadowing,  rather  than  to  the  decreasing  incident  power.  To  investigate  this 
possibility,  we  examined  data  taken  on  the  Ocean  Researcher  1  in  the  South  China  Sea  in  2007  at  a 
wind  speed  of  1 1.7  m/s  and  very  low  grazing  angles.  The  pitching  of  the  ship  varied  the  boresite  of  the 
antennas.  Because  pencil-beam  antennas  were  used,  this  pitching  had  a  large  effect  on  the  signal-to- 
noise  ratio.  In  order  to  determine  the  boresite  during  data  collection,  we  smoothed  the  received  power 
in  range  and  found  the  range  at  which  it  maximized.  Modeling  the  backscatter  using  the  radar  equation 
and  estimates  of  the  NRCS  allowed  us  to  determine  the  antenna  boresite  from  the  range  at  maximum 
received  power.  Figure  4  shows  the  result  of  binning  the  backscatter  into  two  boresite  ranges  that  were 
1°  apart.  In  the  boresite  range  0°  to  0.5°  (Figure  4a),  the  average  SNR  was  higher  and  the  concealed 
fraction  was  lower  than  in  the  boresite  range  1°  to  1.5°  (Figure  4b)  for  every  range  above  1000  m.  Yet 
the  grazing  angle  at  any  particular  range  was  the  same  for  both  boresites.  This  shows  that  the  increase 
in  concealed  fraction  is  due  to  the  decreased  incident  power  for  larger  boresites  and  not  due  to  a 
change  in  grazing  angle.  This  is  true  both  for  HH  and  W  polarization,  again  indicating  partial 
shadowing,  not  geometric  shadowing. 


Figure  3.  Average  SNR  in  dB  divided  by  50  (solid  curve)  and  concealed  fractions  (asterisks)  as 
function  of  range  for  data  taken  on  the  R/V  Thompson  at  W  polarization.  The  dashed  curves  are 
the  predictions  of  Milder’s  geometrical  shadowing  theory  discussed  in  the  text,  a)  Wind  7.3  m/s  from 
9°T,  antenna  boresite  1°  to  2°  looking  toward  7.5°T.  b)  Wind  16.9  m/s  from  341°T,  antenna  boresite 
1°  to  2°  looking  toward  342.5°  T.  Range  resolution  in  7.5  m 
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Figure  4.  Average  SNR  in  dB  divided  by  100  and  concealed  fractions  for  two  different  boresite 
angles  (top  and  bottom)  for  data  taken  on  the  0R1  in  2007.  Wind:  11. 7  m/s  from  25°  T.  a)  Boresite 
=  0.0°  to  0.5°  looking  toward  51°T.  b)  Boresite  =  1.0°  to  1.5°  looking  toward  51°T. 

Range  resolution  is  30  m. 

Figure  5  shows  the  magnitude  and  phase  of  the  coherence  functions  for  HH  and  VV  polarization 
looking  upwind  and  downwind  for  data  taken  in  the  South  China  Sea.  This  data  set  was  collected  with 
a  30  m  range  resolution  so  only  modulation  by  waves  longer  than  60  m  can  be  observed  in  the 
coherence  function.  We  consider  up/down  wave  to  be  up/down  wind.  In  the  figure,  positive  phases 
indicate  that  maxima  of  received  power  modulations  occur  at  shorter  ranges  than  maxima  of  scatterer 
velocity  modulations.  The  figure  shows  a  clear  dependence  of  the  phase  between  received  power  and 
scatterer  velocity  on  both  look  direction  and  polarization. 
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Figure  5.  Magnitude  and  phase  of  the  coherence  function  between  wave-induced  variations  in 
received  power  and  scatterer  velocity.  Magnitudes  are  solid  lines  and  phases  are  asterisks  with  95% 
intervals  shown.  Upwave  and  downwave  measurements  were  made  within  three  hours  of  each 
other.  The  wind  speed  was  10.6  m/s  from  45  °T  and  the  antenna  looked  toward  45  °T  during  the 
upwave  measurements.  The  wind  had  not  changed  but  the  antenna  looked  toward  225  °T  during  the 

downwave  measurements. 
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We  can  model  the  expected  result  of  wave  modulation  on  a  bound  wave/free  wave  picture  similar  to 
that  invoked  earlier  in  both  open  ocean  and  wave  tank  studies  (Plant,  1997;  Plant  et  ah,  1999a,  1999b, 
2004).  The  idea  is  that  long  waves  on  the  ocean  both  modulate  free,  wind-generated  capillaries,  as  the 
standard  composite  surface  model  postulates,  and  cause  intermediate  scale  waves  to  break,  producing 
bound  or  breaking  waves.  The  resultant  backscatter  is  the  sum  of  these  processes.  Both  processes 
vary  with  position  on  the  long  wave  and  we  simplify  the  problem  by  considering  these  variations  only 
to  first  order  in  long  wave  slope.  Then  for  a  sinusoidal  long  ocean  wave, 

r|  =  A  exp{-i(k*R-(ot)} 

we  write  the  (detrended)  wave-modulated  power  and  velocities  as 

Pr  =  <Pr>{Ff  [mh  +  mt]  +  Fbmb}kr| 

Vs=  {Ff+Fbmv}  K-Ckr| 

where  C  is  the  long  wave  phase  speed,  k  is  a  unit  vector  in  the  antenna-look  direction,  <  >  indicates  an 
average  over  the  time  series,  Ff  is  the  fraction  of  backscatter  due  to  free  waves,  Fb  is  the  fraction  of 
backscatter  due  to  bound/breaking  waves,  and  k  =  27t/L  where  L  is  long-wave  length.  Note  that  Ff  +  Fb 
=  1 .  The  different  types  of  modulation  transfer  functions  are  mh,  the  hydrodynamic  MTF  of  the  free 
waves;  mt,  the  tilt  MTF  of  the  free  waves;  mb,  the  bound/breaking  wave  MTF.  The  free  waves  are 
advected  by  the  long  wave  orbital  velocity,  which  is  written  Ckp.  The  transfer  function  mv  accounts 
for  the  fact  that  the  amplitude  and  phase  of  the  bound  wave  velocity  variation  is  not  the  same  as  the 
long  wave  orbital  velocity.  The  amplitude  of  the  bound  wave  velocity  is  the  phase  speed  of  the  parent 
wave  producing  the  bound  waves.  Phases  of  the  modulations  are  with  respect  to  the  long  wave  crest 
and  are  positive  in  the  direction  of  wave  travel. 

The  hydrodynamic  MTF  includes  both  the  modulation  of  the  freely  propagating  short  wave  amplitudes 
and  range  changes  caused  by  the  long  wave  amplitudes;  both  are  in  phase  with  the  long  wave.  The 
orbital  velocity  is  also  in  phase  with  the  long  wave.  The  tilt  mtf  is  always  90°  out  of  phase  with  the 
long  wave  and  occurs  on  the  side  of  the  crest  toward  which  the  antenna  is  looking.  Finally,  Plant 
(1997)  found  that  the  maximum  of  mb  occurred  approximately  130°  in  front  of  the  long  wave  crest.  We 
therefore  model  these  MTFs  as  follows: 

mh  =  2+3/(hk) 

mt  =  15  exp[+br/2],  upwave 

mt  =  15  exp[-iji/2],  downwave 

mb  =  2  exp[i27i(130/180)] 

C=  Vg/k 

mv  =  (1.5/C)  exp[i2;i(130/180)] 
where  h  is  the  height  of  the  antenna. 
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HH  -  Upwave  -  4>p  -  4>v  =  -2°  W  -  Upwave  -  4»p  -  <frv  =  48° 


Figure  6.  Simulated  modulations  of  received  power  and  scatterer  velocity  by  a  sinusoidal  long  wave 
of  1  m  amplitude  and  100  m  wavelength.  The  circles  and  asterisks  are  located  on  the  surface, 
circles  indicating  wind  waves  and  asterisks  indicating  breaking  waves.  The  relative  size  of  the 
circles  and  asterisks  indicates  their  importance  in  HH  and  VV  backscatter.  The  dashed  curve  shows 
the  modulated  velocity,  and  the  dash-dotted  curve  shows  the  modulated  power.  Waves  in  the  upper 
panels  move  to  the  left,  toward  the  radar,  while  those  in  the  lower  panels  move  to  the  right,  away 
from  it.  The  phase  differences  between  power  and  velocity,  <pp  -  <pY,  shown  in  the  figure  agree  with 
the  measurements  shown  in  Figure  10  to  within  about  10°. 


We  may  now  determine  the  behavior  of  the  phase  of  the  correlation  function  for  different  polarizations 
and  look  directions.  We  take  Ff  =  0.9  looking  upwind  and  Ff  =  1  looking  downwind.  For  HH 
polarization  the  bound/breaking  waves  are  dominant  but  a  small  feature  at  free  wave  frequencies  can 
be  seen.  Therefore,  we  take  Ff  =  0.15  looking  upwind  and  Ff  =  0.3  looking  downwind.  We  illustrate 
the  modulation  process  using  a  long  wave  with  A  =  1  m  and  L  =  100  m. 

Figure  6  shows  the  result  of  this  simple,  linear  simulation.  Even  though  the  model  is  surely  too  simple, 
especially  for  the  bound/breaking  waves,  the  phase  differences  between  received  power  and  scatterer 
velocity  modulations  are  close  to  those  shown  in  Figure  10  in  all  cases.  One  can  see  that  power  and 
velocities  maximize  in  most  cases  in  regions  where  geometric  shadowing  would  predict  no  backscatter 
(see  Figure  1).  For  VV  polarization,  the  bound/breaking  waves  make  only  a  minor  difference,  the  shift 
from  ±90°  phases  upwind  and  downwind  being  mostly  due  to  the  hydrodynamic  modulation  of  the  free 
waves.  For  HH  polarization,  on  the  other  hand,  bound/breaking  waves  are  crucial. 
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IMPACT/APPLICATION 

It  is  common  to  think  of  low-grazing-angle  backscatter  from  the  ocean  in  terms  of  geometric 
shadowing.  This  work  shows  that  on  the  open  ocean  this  is  a  very  bad  model.  In  fact,  our  work 
indicates  that  on  the  open  ocean  no  evidence  exists  to  indicate  that  any  reduction  at  all  in  the  incident 
fields  takes  place  in  geometrically  shadowed  regions  at  either  VV  or  HH  polarization.  Therefore 
models  of  low-grazing-angle  backscatter  will  have  to  focus  on  more  realistic  causes  in  regions  of  very 
low  backscatter,  such  as  very  low  scatterer  intensities.  They  cannot  simply  be  brushed  aside  as 
“shadowed  regions’. 

TRANSITIONS 

The  results  of  this  project  have  not  yet  been  transitioned  for  operational  use. 

RELATED  PROJECTS 

The  work  carried  out  in  this  project  to  date  is  closely  related  to  projects  aimed  at  measuring  waves 
around  ships  because  it  addresses  the  causes  of  wave  modulation  at  both  HH  and  VV  polarization.  The 
relevant  projects  are  ONR’s  Ship  Guidance  MURI  and  Hi  Res  DRI  and  the  ongoing,  more  applied 
projects  in  the  Environment  and  Ship  Motion  Forecasting  program. 
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